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ABSTRACT 

Recent observations with the Chandra X-ray Observatory have provided us 
with the capability to discriminate point sources, such as the supermassive black 
hole Sgr A*, from the diffuse emission within the inner 10" of the Galaxy. The 
hot plasma producing the diffuse X-radiation, estimated at ~ 7.6 x 10 31 ergs s _1 
arcsec -2 in the 2-10 keV band, has a RMS electron density rs 26 cm -3 and a 
temperature kT w 1.3 keV, with a total inferred mass of ~ 0.1 M & . At least some 
of this gas must be injected into the ISM via stellar winds. In the most recent 
census, about 25 bright, young stars have been identified as the dominant sources 
of the overall mass efflux from the Galactic center. In this paper, we use detailed 
3-dimensional SPH simulations to study the wind-wind interactions occurring in 
the inner 3 parsecs of the Galaxy, with a goal of understanding what fraction, if 
any, of the diffuse X-ray flux measured by Chandra results from the ensuing shock 
heating of the ambient medium. We conclude that this process alone can account 
for the entire X-ray flux observed by Chandra in the inner 10" of the Galaxy. 
Understanding the X-ray morphology of the environment surrounding Sgr A* 
will ultimately provide us with a greater precision in modeling the accretion of 
gas onto this object, which appears to be relatively underluminous compared to 
its brethren in the nuclei of other galaxies. 

Subject headings: accretion — black hole physics — Galaxy: center — radiation mech- 
anisms: thermal — stars: winds — X-rays: diffuse 

1. Introduction 

The inner few parsecs of the Galaxy contain several observed components that coexist 
within the central deep gravitational potential well. The pointlike, nonthermal radio source 
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Sgr A* appears to be the radiative manifestation of a ~ 3 million M concentration of dark 
matter within ~ 0.015 pc of the nucleus. This concentration, probably a central black hole, 
is surrounded by a cluster of evolved and young stars, a molecular dusty ring, ionized gas 
streamers, diffuse hot gas and Sgr A East, a shell-like, nonthermal radio source believed to 
be the remnant (SNR) of a supernova explosion at the center of the Galaxy some 10, 000 
years ago (see Melia & Falcke 2001). 

Observations at infrared and, particularly, radio wavelengths have been very effective 
in unraveling the complex behavior of the central components through their mutual interac- 
tions. But only through recent observations with the Chandra X-ray Observatory have we 
been able to study the X-ray emission from Sgr A* and its surroundings with sub-arcsecond 
resolution and a broad X-ray band imaging detector. These capabilities have provided us 
with the means to discriminate point sources, such as Sgr A*, from the diffuse emission due 
to hot plasma in the ambient medium. 

Baganoff et al. (2003) report that the Chandra/ AClS-l 0.5-7 keV image of the inner 
17' x 17' region of the galaxy contains remarkable structure; its detail is sufficient to allow 
comparisons with features seen in the radio and IR wavebands. The diffuse X-ray emission is 
strongest in the center of Sgr A East. Based on these X-ray observations, Maeda et al. (2001) 
conclude that Sgr A East is a rare type of metal-rich, "mixed-morphology" SNR that was 
produced by the Type II explosion of a 13-20 M progenitor. The X-ray emission from Sgr 
A East appears to be concentrated within the central 2-3 pc of the 6 x 9 pc radio shell and 
is offset by about 2 pc from Sgr A* itself. The thermal plasma at the center of Sgr A East 
has a temperature kT 2 keV, and appears to have strongly enhanced metal abundances 
with elemental stratification. 

However, this is not the sole component of hot gas surrounding Sgr A*. Chandra 
has also detected X-ray emission extending perpendicular to the Galactic plane in both 
directions through the position of Sgr A*, which would seem to indicate the presence of a 
hot, "bipolar" outflow from the nucleus. And an additional component of X-ray emitting gas 
appears to have been found within the molecular ring, much closer to Sgr A*. Indeed, the 
western boundary of the brightest diffuse X-ray emission within 2 pc of the Galactic center 
coincides very tightly with the shape of the Western Arc of the thermal radio source Sgr A 
West. Given that the Western Arc is believed to be the ionized inner edge of the molecular 
ring, this coincident morphology suggests that the brightest X-ray emitting plasma is being 
confined by the western side of the gaseous torus. In contrast, the emission continues rather 
smoothly into the heart of Sgr A East toward the eastern side. 

A detailed fit of the emission from the hot gas within 10" of Sgr A* shows that it is clearly 
not an extension of the X-ray emitting plasma at the center of Sgr A East. The inferred 
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2-10 keV flux and luminosity of the local diffuse emission are (1.9 ±0.1) x 10~ 15 ergs cm -2 
s _1 arcsec -2 and (7.6+ 2 ;g) x 10 31 ergs s" 1 arcsec -2 , respectively. Based on the parameters 
of the best-fit model, it is estimated that the local, hot diffuse plasma has a RMS electron 
density (h 2 ) 1 / 2 ~ 26 cm -3 . Assuming that this plasma has a filling factor of unity and is 
fully ionized with twice solar abundances, its total inferred mass is ^i oca i ~ 0.1 M . The 
local plasma around Sgr A* appears to be distinctly cooler than the Sgr A East gas, with a 
temperature kT = 1.3+q; 2 keV. By comparison, the best-fit model for Sgr A East indicates 
that the plasma there has about 4 times solar abundances, and twice this temperature (see 
above). 

At least some of this local, hot plasma must be injected into the ISM via stellar winds. 
There is ample observational evidence in this region for the existence of rather strong outflows 
in and around the nucleus. Measurements of high velocities associated with IR sources in 
Sgr A West (Krabbe et al. 1991) and in IRS 16 (Geballe et al. 1991), the H 2 emission in the 
molecular ring and from molecular gas being shocked by a nuclear mass outflow (Gatley et 
al. 1986; but see Jackson et al. 1993 for the potential importance of UV photodissociation 
in promoting this H 2 emission), broad Bra, Br7 and He I emission lines from the vicinity 
of IRS 16 (Hall, Kleinmann & Scoville 1982; Allen, Hyland & Hillier 1990; Geballe et al. 
1991), and radio continuum observations of IRS 7 (Yusef-Zadeh & Melia 1991), provide clear 
evidence of a hypersonic wind, with a velocity v w ~ 500-1,000 km s _1 , a number density 
n w ~ 10 3 ~ 4 cm -3 , and a total mass loss rate M w ~ 3-4 xlO -3 M yr -1 , pervading the inner 
parsec of the Galaxy. 

Chevalier (1992) discussed the properties of stellar winds in the Galactic center and 
their ability to produce X-rays, under the assumption that the colliding winds create a 
uniform medium with a correspondingly uniform mass and power input within the central 
0.8 pc of the Galaxy. While the assumption of uniformity may be an over-simplification, 
steady-state models yield central conditions similar to those calculated in our simulations 
and (apparently) observed in nature. 

The possibility that wind-wind collisions in the dense stellar core of the Galaxy might 
produce an observable X-ray signature was considered more recently in the context of tran- 
sient phenomena by Ozernoy et al. (1997). These authors argued that the strong winds from 
OB stars or WR stars would produce shocks that heat the gas up to a temperature T ~ 10 7 
K. Instabilities in the overall outflow could therefore produce transient high temperature 
regions, in which the X-ray luminosity could be as high as 10 33 -10 35 ergs s -1 . 

A preliminary, low-resolution 3-dimensional hydro dynamical simulation by Melia & 
Coker (1999) subsequently revealed that although the tesselated pattern of wind-wind col- 
lision shocks can shift on a local dynamical time scale, the overall concentration of shocks 
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in the central region remains rather steady. This is hardly surprising, given that 25 or so 
dominant wind sources are distributed rather evenly about the center. 

In this paper, we study in detail the emission characteristics of shocked gas produced 
in wind-wind collisions within 3 parsecs of Sgr A*. We present the results of comprehensive, 
high-resolution numerical simulations of the wind-wind interactions, using the latest suite 
of stellar wind sources with their currently inferred wind velocities and outflow rates. We 
use their projected positions and vary the radial coordinates to determine the effect of 
this uncertainty on the X-ray emissivity. In addition, it is clear now from the Chandra 
image that the stellar winds are also interacting directly with the molecular ring, and so we 
incorporate the latter into the set of boundary conditions for our simulations. The predicted 
X-ray luminosity and spatial profile may then be compared with the Chandra data. We 
describe our numerical technique, including the description of this molecular ring and the 
characteristics of the wind sources in § 2 of the paper. In § 3 we present our results, and 
discuss their implications for the conditions in the Galactic center in § 4. 



2. The Physical Setup 

Our calculations use the 3-dimensional smoothed particle hydrodynamics (SPH) code 
discussed in Fryer & Warren (2002), and Warren et al. (2003). The Lagrangean nature 
of SPH allows us to concentrate spatial resolution near shocks and capture the complex 
time-dependent structure of the shocks in the Galactic center. 

We assume that the gas behaves as an ideal gas, using a gamma-law (7 = 5/3) equation 
of state. Since the mass of gas in the simulation at any time is many orders of magnitude 
less than the mass of the black hole and its halo of stars, we calculate gravitational effects 
using only the gravitational potential of the assumed central 2.6 x 10 6 M point mass. 

Each of our two calculations uses slightly different initial conditions; details of each 
calculation are summarized in Table 1. The initial conditions for each simulation assume 
that the Galactic center has been cleared of all mass (perhaps by the supernova explosion 
that produced Sgr A East). Massive stars then inject mass into the volume of solution via 
winds, gradually increasing the mass of gas in the Galactic center. The number of particles 
in each simulation initially grows rapidly, but reaches a steady state (~ 7 million particles) 
when the addition of particles from wind sources is compensated by the particle loss as 
particles flow out of the computational domain or onto the central black hole. The particle 
masses vary from 9.4 x 10" 7 M to 1.0 x 10~ 5 M . 

In addition, we place a torus of molecular material in a circumnuclear disk surrounding 
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the inner 1.2 pc of the Galaxy. The boundary conditions (outer and central black hole), 
the circumnuclear disk, and the stellar wind sources are the only modifications made to 
the basic SPH code of Warren et al. (2003) to run these simulations, and we discuss these 
modifications in more detail in the next three subsections. 



2.1. Boundary Conditions 

The boundary conditions are the major particle sinks in the simulation. We assume 
that once these particles have achieved escape velocity beyond a certain distance, nothing 
confines the material flowing out of the Galactic center. The computational domain is a cube 
approximately 6 pc on a side, centered on the black hole. To simulate "flow-out" conditions, 
particles passing through this outer boundary are removed from the calculation. 

In addition, particles passing through an inner spherical boundary with a radius of 
1.9 x 10 17 cm (equivalent to 2.5 x 10 5 R s or ~ 0.7 Bondi-Hoyle capture radii for Sgr A*) 
centered on the black hole are also removed. This assumes that the accretion onto Sgr A* 
continues roughly at free-fall below this point. Such a simplified assumption only provides a 
limit on the X-ray luminosity of the point-source of Sgr A*, but it does not affect the diffuse 
X-ray emission in the Galactic center that is the central topic of this paper. 

2.2. Circumnuclear Disk 

HCN synthesis data, taken at 5" — 10" spatial and 4 km s _1 spectral resolution, reveal 
a highly inclined, clumpy ring of molecular gas surrounding the ionized central ~ 2 pc of 
the Galaxy (Giisten et al. 1987). This ring seems to be the inner edge of a thin molecular 
structure extending out to ~ 5 — 7 pc from the center. Within this cavity, and orbiting about 
Sgr A*, is the huge H II region known as Sgr A West comprised of a 3-armed mini-spiral 
(Lo & Claussen 1983; Lacy et al. 1991). The molecular and ionized gas in the central cavity 
appear to be coupled, at least along the so-called western arc (one of the spiral arms), which 
appears to be the ionized inner surface of the molecular ring. It has been suggested that the 
northern and eastern arms may be streamers of ionized gas falling toward Sgr A* from other 
portions of the molecular ring (see Melia & Falcke 2001, and references cited therein). 

Since 1987, observations of various tracers of molecular gas (notably H 2 , CO, CS, O 
I, HCN, NH 3 ) have somewhat refined our view of the circumnuclear disk (or CND, as it 
is sometimes called), leaving us with the following set of defining characteristics: a total 
mass of ~ 3 x 10 4 M Q , a very clumpy distribution with an estimated volume filling factor 
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of ~ 0.01, a typical clump mass of ~ 30 M Q , size ~ 0.1 pc, and temperature ~ 100 K (see 
Gatley et al. 1986; Serabyn et al. 1986; Giisten et al. 1987; DePoy et al. 1989; Sutton et al. 
1990; Jackson et al. 1993; Marr et al. 1993). As many as 500 clumps fill its overall structure. 
Its height at the inner edge is ~ 0.5 pc, and it flares outward with distance from the center. 

The Chandra 0.5 — 7 keV image of the central 1.3' x 1.5' region of the Galaxy (see Fig. 4 
in Baganoff et al. 2003), overlaid on the VLA 6-cm contours of Sgr A* and Sgr A West, shows 
that the western boundary of the brightest diffuse X-ray emission coincides very closely with 
the shape of the western arc, adding credence to the notion that the latter tracks the inner 
edge of the CND if one adopts the view that the X-ray glowing plasma is itself being confined 
by the western side of the molecular ring. 

The CND therefore appears to be an essential element in the (partial) confinement 
of the hot gas, and we incorporate it into our simulations using a simple approach that 
nonetheless retains the CND's important features. To keep the volume of solution tractable, 
we do not model the entire ~ 14-pc disk, but rather introduce an azimuthally-symmetric 
torus with inner radius 1.2 pc (see Giisten et al. 1987) and thickness 1 pc. The latter is a 
fair representation of the CND's observed height at this radius. This structure does not, of 
course, account for the CND's outer regions, but as we shall see, the hot, X-ray emitting gas 
does not penetrate far into the molecular boundary, so the outer CND is not directly relevant 
to this simulation. We assume that this torus has a mass of ~ 10 4 M & , and is composed of 
200 spherical clumps, each of mass 50 M & . We also incline it by 60° from the plane of the 
sky, orienting its principal axis (in projection) along the Galactic plane. 

2.3. Wind Sources 

It is well known by now that the Galactic center wind is unlikely to be uniform, since 
many stars contribute to the mass ejection (see, e.g., Melia & Coker 1999). For these 
calculations, we assume that the early-type stars enclosed (in projection) within the Western 
Arc, the Northern Arm, and the Bar produce the observed wind. Thus far, 25 such stars 
have been identified (Genzel et al. 1996), though the stellar wind characteristics of only 8 
have been determined from their He I line emission (Najarro et al. 1997). Figure 1 shows 
the positions (relative to Sgr A*; 1" 0.04 pc at the Galactic center) of these wind sources; 
the size of the circle marking each position corresponds to the relative mass loss rate (on a 
linear scale) for that star. IRS 13E1 and IRS 7W seem to dominate the mass outflow with 
their high wind velocity (~ 1,000 km s -1 ) and a mass loss rate of more than 2 x 10" 4 M 
yr _1 each. 
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Wind sources are modeled as literal sources of SPH particles. New particles are added in 
shells around each wind source as the simulation progresses and the existing particles move 
outward from each source. Using wind velocities inferred from observations (Najarro et al. 
1997), shells around each wind source are added at a rate such that the spacing between 
shells is approximately equal to the spacing between particles within each incoming shell. 
Masses of the incoming particles are chosen to match the inferred mass loss rate for each 
source. The stars without any observed He I line emission are assigned a wind velocity of 
750 km s _1 and an equal mass loss rate chosen such that the total mass ejected by the stars 
used here is equal to 3 x 1CT 3 M yr _1 , the overall mass outflow inferred for the Galactic 
center region (see, e.g. Melia 1994). The initial temperature of the stellar winds is not 
well known; for simplicity, we assume here that all the winds are Mach 30 at their point of 
ejection from the stellar surface. This corresponds to a gas temperature of ~ 10 4 K. We note, 
however, that our results are insensitive to this value, since the temperature of the shocked 
gas is set primarily by the kinetic energy flux carried into the collision by the winds. For 
the calculations reported here, the sources are assumed to be stationary over the duration 
of the simulation. 

An additional uncertainty is the location in z (i.e., along the line of sight) of the wind 
sources. To test the sensitivity of our results to the choice of z, we carry out two simulations 
with different coordinate assignments, and compare the overall shocked gas configurations 
and their X-ray emissivities. For the first choice of z positions — listed as z x in Table 2 and 
used in Simulation 1 — we determine this third spatial coordinate randomly, subject to the 
condition that the overall distribution in this direction matches that in x and y. With this 
proviso, all these early-type stars are located within the central few parsecs surrounding Sgr 
A*. For the second choice of positions — listed as Z2 and used in Simulation 2 — we attempt 
to calculate a reasonable upper limit on the luminosity produced in the central parsec by 
moving all the sources closer together in z, toward the z = plane. 



3. The Gas Profile and Spectrum 

Wind-wind collisions in the central parsec create a complex configuration of shocks that 
efficiently converts the kinetic energy of the winds into internal energy of the gas. Figure 2 
shows isosurfaces of specific internal energy in the central cubic parsec of Simulation 1, 
~ 10, 000 years after the beginning of the calculation. The blue surfaces indicate regions of 
gas with low specific internal energy, which lie near the wind sources; the red surfaces mark 
regions of high specific internal energy, where gas has passed through multiple shocks. 26% 
of the total energy in the central parsec has been converted to internal energy via multiple 
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shocks: the total kinetic energy of material in the central parsec is 7.7 x 10 48 ergs, while the 
total internal energy is 2.7 x 10 48 ergs. 

In order to calculate the observed continuum spectrum, we assume that the observer 
is positioned along the negative z-axis at infinity and we sum the emission from all of 
the winds — and their shocks — produced by the 25 stars in Table 2. For the conditions we 
encounter here, scattering is negligible and the optical depth is always less than unity. For 
these temperatures and densities, the dominant components of the continuum emissivity are 
electron-ion (e e j) and electron-electron (e ee ) bremsstrahlung. 

Figures 3 and 4 show the 2-10 keV luminosity per arcsec 2 from the central 10" of 
Simulations 1 and 2, respectively, versus time since the beginning of each calculation. The 
winds fill the volume of solution after ~ 4000 years, which coincides with the point in each 
figure where the luminosity stops rising rapidly and begins varying by only ~ 5% around a 
steady value. The inset plot in each figure shows variation of the luminosity on a timescale 
of ~ 10 years. Data points in the insets are plotted for each of 990 time steps, where the size 
of each time step is 0.16 years. Variation between consecutive data points is due primarily 
to numerical noise, but variations on timescales of several years indicate noticeable changes 
in the temperature and density of the gas in the central 10". 

The range of values of the average luminosity per arcsec 2 measured in the central 10" 
of our simulations illustrates the sensitivity of the luminosity to the positions of the wind 
sources. Luminosities for each of our simulations are reported in Table 1; for comparison, 
Baganoff et al. reported a 2-10 keV luminosity of (7.6l 2 ;g) x 10 31 ergs s" 1 arcsec -2 from the 
local diffuse emission within 10" of Sgr A*. Results from both simulations fall within the 
error bars reported by Baganoff et al. 

Table 1 also includes two different estimates for the temperature of the plasma within 
10" of Sgr A*. T m ass is an average temperature calculated using the mass of each particle 
as a weighting factor, while Ti um uses the luminosity of each particle as a weighting factor. 
Since the temperature is determined primarily by the kinetic energy flux initially carried by 
the winds, and since the initial wind velocities in each calculation are the same, the average 
temperature for the plasma at the center of each calculation is essentially the same. 

Figures 5 and 6 show the average luminosity per volume and average luminosity-weighted 
temperature as a function of radius from the center of Simulation 1. The majority (77% of 
the overall 2-10 keV luminosity) of the 2-10 keV emission comes from the central 0.4 pc, 
where it is produced in wind-wind collisions. There is also noticeable emission (2.5% of the 
overall luminosity) from a high-temperature region centered on 1.2 pc, which is the location 
of the inner edge of the CND; emission here is produced when winds collide with the clumps 
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in our model torus. 

Figures 7 and 8 show contours of 2-10 keV luminosity per arcsec 2 from Simulations 1 
and 2, respectively, 10, 000 years after the winds were turned on. The dotted circle overlaid 
on each figure shows the size of the 10" central region. The zones of relatively high luminosity 
at the north-east and south-west corners of each figure, and along the northwest side of the 
Simulation 2 contours, result from collisions of the outflowing wind with the inner edge of 
the CND. 

The gas distribution in a multiple-wind source environment like that modeled here is 
distinctly different from that of a uniform flow past a central accretor (Coker & Melia 1997). 
Unlike the latter, the former does not produce a large-scale bow shock, and therefore the 
environmental impact of the gravitational focusing by the central dark mass has significantly 
less order in this case. We shall defer a more extensive discussion of this point to a later 
publication, in which we report the results of a multiple-wind source simulation for the 
accretion of shocked gas onto the central, massive black hole. 

4. Discussion 

Using only the latest observed stellar positions and inferred mass loss rates and wind 
velocities, we have produced self-consistent X-ray luminosity maps of the Galactic center. 
It appears that the diffuse X-ray luminosity and temperature within 10" of Sgr A* can be 
explained entirely by shocked winds; the luminosities from both of our simulations fall within 
the error bars associated with the diffuse 2-10 keV X-ray luminosity measured by Chandra. 
This seems rather significant in view of the fact that, were the diffuse X-rays surrounding 
Sgr A* being produced by another mechanism, we would need to lower the overall stellar 
mass-loss rate or wind velocity by a factor of 2 or more in order to render the emissivity 
produced by wind-wind collisions unobservable. The luminosities that we calculate are fairly 
insensitive to the placement of the wind sources along the line of sight to the Galactic Center; 
the 2-10 keV luminosity per arcsec 2 differs by only ~ 15% between Simulation 1 — in which 
sources are randomly distributed in z — and Simulation 2 — in which sources are compressed 
toward Sgr A* perpendicular to the plane of the sky. 

Although the wind mass-loss rate is 3 x 10 -3 M & yr _1 , the average accretion rate though 
our inner boundary is only 4.25 x 10 -4 M yr _1 for Simulation 1 and 4.59 x 10~ 4 M yr _1 
for Simulation 2. Even so, this is much higher than what one would expect from the X-ray 
emission at this inner boundary. Since we do not model this accretion, we have neglected 
a number of effects from magnetic fields to radiation pressure. We postpone this detailed 
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analysis to a future paper. 

One possible explanation for the relatively low luminosity of Sgr A* has been that the 
supernova that created Sgr A East swept most of the gas out of the environment surrounding 
Sgr A*. However, our calculations show that the stellar winds within 3 parsecs of the Galactic 
center bring the environment near the black hole back to a steady state within ~ 4000 years. 
In addition, the gas temperature and X-ray luminosity calculated from our results add weight 
to Baganoff et al.'s hypothesis that the X-ray emitting gas surrounding Sgr A* really is a 
different component than that at the center of Sgr A East, which has a temperature closer 
to 2 keV. 

Collision of the wind with the inner edge of the CND may produce the western arm 
of the minispiral. Our simple model for the CND produces clear emission at the location 
of the inner edge of the molecular ring, although we do not reproduce the entire complex 
morphology of the X-ray emitting gas seen by Chandra. 

The excellent agreement between our calculated X-ray luminosities and the value of the 
luminosity measured by Chandra solidifies our understanding of the gas dynamics in the 
Galactic center. Future calculations will move the inner boundary of the simulation inward 
and explore the dynamics of the shocked gas as it falls toward the central black hole. 
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Table 1. Simulation Properties 



Simulation 


z positions 


10" L 


T 

J - mass 








(ergs s _1 arcsec -2 ) 


(keV) 


(keV) 


SI 


Zl 


6.45 x 10 31 


1.72 


0.25 


S2 


Z2 


7.50 x 10 31 


1.76 


0.25 
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Table 2. Parameters for Galactic Center Wind Sources 



Star 


x a 




zi a 




V 


M 


(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


/l 1 \ 

(km s J 


(10- 5 M yr- 1 ) 
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Fig. 1. — Location of the 25 wind-producing stars used in the simulations reported here, 
relative to the position of Sgr A* indicated by the * symbol. The radius of each circle 
corresponds (on a linear scale) to that star's mass loss rate. Setting the scale is 13E1, with 
M = 7.9 x 1CT 4 M yr" 1 . 
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Fig. 2. — Isosurfaces of specific internal energy from Simulation 1, 10, 000 years after the 
winds were turned on. The blue surfaces correspond to a specific internal energy of 2.5 x 
10 12 ergs g _1 ; the red surfaces correspond to 3.8 x 10 15 ergs g -1 . 
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Fig. 3. — 2-10 keV X-ray luminosity per arcsec 2 from the central 10" versus time for Simu- 
lation 1. The large plot shows the variation in luminosity over the entire calculation, while 
the inset plot shows variation on a timescale of ~ 10 years. The winds fill the volume of 
solution after ~ 4000 years. The dashed line indicates the luminosity measured by Chandra, 
while the dotted line indicates the value of the lower error bar on that measurement. The 
upper error bar, at 1.02 x 10 32 ergs s -1 arcsec -2 , is not visible on this graph. 



-18- 



8xl0 31 - 



6xl0 31 



4xl0 31 - 



2xl0 31 - 



- 




7.6xl0 31 - 



7.4xl0 31 11111 



j i i i i_ 



9600 9650 



9700 



_L 



2000 



4000 



6000 



8000 



10000 



t (yr) 



Fig. 4. — 2-10 keV X-ray luminosity per arcsec 2 from the central 10" versus time for Simu- 
lation 2. The large plot shows the variation in luminosity over the entire calculation, while 
the inset plot shows variation on a timescale of ~ 10 years. The winds fill the volume of 
solution after ~ 4000 years. The dashed line indicates the luminosity measured by Chandra, 
while the dotted line indicates the value of the lower error bar on that measurement. The 
upper error bar, at 1.02 x 10 32 ergs s -1 arcsec -2 , is not visible on this graph. 
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Fig. 5. — 2-10 keV X-ray luminosity per volume per 0.035 pc radial "bin" versus distance 
from the central black hole in Simulation 1. Collisions between the winds and the inner edge 
of the CND are visible as a bump near 1.2 pc. 
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Fig. 6. — Average luminosity-weighted temperature per 0.035 pc radial "bin" versus distance 
from the central black hole in Simulation 1. Collisions between the winds and the inner edge 
of the CND are visible as a spike near 1 pc. 
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Fig. 7. — Contours of column-integrated 2-10 keV X-ray luminosity per "bin" from Simu- 
lation 1. There are 100 bins along each axis; each bin is 0.9" by 0.9". In order from black 
to red to green to blue, the luminosities indicated by the contours are 10 32 ergs s _1 per bin, 
10 31 ergs s^ 1 per bin, 10 30 ergs s -1 per bin, and 1.5 x 10 29 ergs s _1 per bin. The dashed 
circle indicates the extent of the "local" region around Sgr A*; it has a radius of 10" and is 
centered on Sgr A*. 
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Fig. 8. — Contours of column-integrated 2-10 keV X-ray luminosity per "bin" from Simu- 
lation 2. There are 100 bins along each axis; each bin is 0.9" by 0.9". In order from black 
to red to green to blue, the luminosities indicated by the contours are 10 32 ergs s _1 per bin, 
10 31 ergs s^ 1 per bin, 10 30 ergs s -1 per bin, and 1.5 x 10 29 ergs s _1 per bin. The dashed 
circle indicates the extent of the "local" region around Sgr A*; it has a radius of 10" and is 
centered on Sgr A*. 



